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Application of Two Dimensional Unsteady 
Aerodynamics to a Free-Tip Rotor Response Analysis 


Leslie Yates 

Russmark Inc. 
602 Stendhal Ln. 
Cupertino, C A 95014 


Hiroyuki Kumagai 

The University of Kansas Center of Research Inc. 
2291 Irving Hill I)r.. Campus West 
Lawrence. KS 66015-2969 


Summary 

The fret-tip rotor utilizes a rotor blade tip which is structurally decoupuled from the 
Made inboard section. The tip is fret- to pitch about its own pitch axis to respond to the 
loca' flow angularity changes. The tip also experiences the heaving motion due to the 
flapping of the rotor blade. For an airfoil in any pitching and heaving motion which can 
In expanded into a Fourier series, the lilt and moment calculated by Theodorsen’s theory 
is -imply the linear combination ol the lift and moment calculated for each harmonic. 
These lit and moment are then used to determine the response of the free-tip rotor. A 
parametric study is performed to determine the effect of mechanical damping, mechanical 
spring, -weep, friction, and a constant control moment on the free-tip rotor response 
characterisin'- and the resulting azimuthal lift distributions. The results showed that the 
Irer-tip ha- the capability to suppress the oscillatory lift distribution around the azimuth 
and to eliminate a significant negative lift peak on the advancing lip. This result agrees 
with the result of the previous analysis based on the steady aerody namics. 







1. Introduction 


For a helicopter in forward flight, there is a significant negative lift on the tip of 
the advancing blade. This combined with a positive lift on the retreating side causes an 
oscillatory lift loading which acts as a source of vibration. It has been suggested that a 
free-lip configuration could suppress the oscillatory azimuthal lift distributions at the tip 
and minimize this source of the vibration. 

A. free-tip rotor consists of a conventional inboard section with the outer five to fifteen 
percent of the blade allowed to pilch independently of the inner section of the blade (Figure 

1) . A more detailed discussion can be found in Reference J and 2. In programs FTR2 and 
HOVER, which compute the tip response in forward flight and the tip response to an air jet 
in hover respectively. Kumagai used a steady calculation for the lift coefficient, a constant 
moment coefficient at the 25% chord, and an aerodynamic damping factor calculated from 
the Quasi Vortex Lattice Method to stud) the behavior of the free-tip (Reference 1 and 

2) . Mis study indicated that the free-tip does indeed make the azimuthal lift distributions 
less oscillatory. 

Hoeing V'orlol also developed a program to determine the behavior of the free-tip. In 
their program. B-CS5. the lift and moment are calculated from Theodorsen’s theory for wings 
in pitching and heaving motion (Reference 3). The aerodynamic damping is automatically 
included in these calculations. In B-K5 it is assumed that the local angle of attack, the 
inflow angle, the flapping angle, and the downwash velocity can all be approximated by a 
constant term and a sinusoidal term with the angular velocity of the rotor. 

The purpose of the present study is to include Theodorsen's theory in the existing 
two dimensional numerical models. FTR2 and HOVER, and investigate its effects on the 
free-tip response characteristics. Although no modifications were made to the lift and 
moment derived from Theodorsen's theory, the local angle of attack, the inflow angle, the 
flapping angle, and the downwash velocity are all expanded into Fourier series, and the lift 
and moment are then uniquely determined for each harmonic. The computer programs 
developed to determine the response of the free-tip were UNSTEADY and rNSTHOVR. 

For a small scale model in forward flight, the local angle of attack was calculated 
using three methods. FTR2, B-G5, and UNSTEADY. The differences in the solutions are 
discussed in reference to the formulation of the three methods. For the same model in 
hover, the response to an air jet was examined using HOVER and IWSTIIOVR. 

For the small scale model and a full scale model, the el feet of changing a number of 
parameters on the azimuthal lift distribution was studied. The parameters varied were the 
constant moment coefficient at the 25% chord, the constant control moment, the sweep in 
terms of the aerodynamic center offset of the free-lip from its pitch axis, the friction, the 
mechanical damping, and the mechanical spring constant. 

•> 




2. List of Symbols 





A On 

A c „m 
Ae„m 
^rt Bn 

b o 
be 
c<i 
c\ 

Cl 

CL 

CL 

Cm 

C„ 



./ 

k, 

l^e 


Fourier coefficients of £(1) 

Fourier coefficients of £(cos nut) 

Fourier coefficients of £(sinnu>/) 

harmonic coefficients of a 

aerodynamic damping coefficient {lb ft sec/ ft) 

mechanical damping coefficient (lb ft sec/ ft) 

drag coefficient 

lift coefficient 

lift curve slope (/rad) 

lift coefficient due to the steady part of the a, 0. -ft . and d> 
lift coefficient due to the n-th harmonic 
moment coefficient 
chord length (ft) 

location of the pitch axis measured from the leading edge in fraction of c (1 

rotor thrust coefficient 

rotor propulsive force coefficient 

down wash velocity (ft/scc) 

harmonic coefficients of the downwash velocity (ft /sec) 
vertical location of the free-lip due to the heaving motion (ft) 
pure imaginary number 

moment of inertia for the free-tip per unit span (slug ft 2 / ft) 
aerodynamic spring coefficient (lb ft/ ft) 
mechanical spring coefficient (lb ft/ ft) 

3 



c 

^■a, (f 

i'ti)' 

TTI 

M 

M 

Mu, Mem, ^ Am 

m,„ m j 
Mdh 
Mj 
M. 

Mj. 

V 

Poe 

9 

9<* 

R 

t 

V 

Vo, 
ti, t', to 

V 

x,y.z 


operator acting upon a, defined by eq. (27) 

reference lengths used in calculating the feathering moment (ft) 

span of the free-tip (ft) 

mass of the free-tip (slug/ ft) 

Mach number 

operator acting upon </>, ^ and [i , defined by eq. (27) 

Fourier coefficients of M 

moment measured at the quarter chord (Id ft/ ft) 

forward and aft masses contributing to the feathering moment (lb ft/ ft) 

moment due to drag and the dihedral (lb ft/ ft) 

moment caused by friction (lb ft/ ft) 

moment caused by the aerodynamic spring (lb ft! ft) 

feathering moment (lb ft/ ft) 

total pressure (lb/ ft 2 ) 

total pressure at infinity ( lb/ ft 2 ) 

dynamic pressure (lb/ ft 2 ) 

constant control moment (lb ft/ ft) 

radius .of the rotor (ft) 

time (sec) 

velocity of the rotor blade tip (fl/scc) 
free stream velocity (fl/scc) 
perturbation velocity (ft j sec) 
forward velocity of the helicopter (ft /sec) 

Cartesian coordinates (ft) 



location of a shed vortex (ft) ' 

location of the pitch axis measured from mid-chord of the airfoil (ft) 
local angle of attack (rad) 

constant term in the Fourier expansion of the local angle of attack (rad) 
flapping angle (rad) 

effective pilch axis offset, distance between pilch axis and effective 
aerodynamic center in fraction of the chord 

strength of shed vortex 

harmonic coefficients of AT 

free-tip dihedral, positive down (rad) 

pitch angle of the free-tip measured from disk plane (rad) 

rotor advance ratio 

friction coefficient 

free stream air density (slug/ ft*) 

rotor solidity 

velocity potential 

perturbation velocity potential 

inflow angle (rad) 

harmonic, coefficients of the inflow angle (rad) 

azimuthal angle (rad) 

angular velocity of the rotor (rad /.sec) 

angular velocity of airfoil pitching motion ( rad/sec ) 



3. Formulations 

3.1. Theodorsen’s Lift and Moment Coefficients 


For an airfoil in potential flow, the velocity potential can be defined as (Figure 2) 

= * + $' (1) 

where VooX is the free stream velocity potential and <i> # is the perturbation velocity poten- 
tial. 

The pressure distribution on the surface of the airfoil is given by the unsteady, incom- 
pressible Bernoulli’s equation, 


P~ Poo 




Pco dt 2 


2 

2 - oo 


( 2 ) 


The vector (u,v,w) is the perturbation velocity vector. If the perturbation velocities are 
small with respect to the free stream velocity, this equation can be linearized to give 


P-Poo d*' ,, 

= ~ ~7T, — v ^ u 

Poo ot 


(3) 


The lifting pressure is the difference in the pressure on the upper and lower surfaces of the 
wing, 

A „ (d< rr d <\ 

Ap=- 2 />»(- 5 r +t/oo 1 j-j m 

where 


— ^ I upper purf;u 


= |lo 


ower piirCiee 


(3) 


The perturbation velocity potential for an airfoil in pitching and heaving motion 
consists of three parts, the velocity potential due to the pitching motion, $(0), the velocity 
potential due to the heaving motion, and the velocity potential due to shed vortices, 
^(AF). Following Theodorsen's derivation, these potentials evaluated on the surface of 
the airfoil are given by 


$(0) =V x 0VT-s-> t ^ ^ \/7- x 


*(h)=f t \/T Z ~x' 

$ (AF) - T ~2n~^ tan 


( 6 ) 


i ( v'i'-rV ar, 2 ,- A 

\ 1 - xx ( , ) 
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where 0 is pitch angle, x is the location on the airfoil, ^ is the downwash velocity, and 
$(Ar) is the velocity potential for a vortex of strength AT at x 0 (Figure 2). 

To meet the Kutta condition at the trailing edge, the perturbation potential must 
satisfy 

= finite (7) 


dx 


*=j 


Letting + 4>(/i) + $(AT) and substituting into Equation (7) gives 

-x (it a dO (x \ dh\ 


x/T 


+ Ar(l0 > : lrr‘' 10 


w 




— finite 


lx=J 


For the above expression to be finite, the coefficient of must be zero at x = 1. This 

defines an integral equation for the vortex strength density Af(xo). 


1 f™ x *0 + 1 rr a dO f 1 \ dh 

sj,. ArM v^‘' II, = t, “" + *( 2 - x v + 5r 




In the helicopter environment it is possible to expand 0 and into the Fourier series 


oo 


0=0 o +^0 n e 1n 


Ut 


rj — 1 
DO 


, V d tjle' nUt 

dl dl *-*• dl 

. * • n=. I • 


(10) 


where fl is the angular velocity of the rotor and the coefficients 0 n and are complex 
to allow for the phase angles. 

Assuming that the vortices are carried away with the free stream velocity, the vortex 
density can be written as 

/ T-s \ 

(H) 


Ar(x o ,0 = Ar (t - ~) 

This also can be expanded into a Fourier series, 


OO 


Ar(*„,r) - Ar 0 + ^2 


/<') 


(12) 


n~ 1 
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It should be noted that the “free stream” velocity for the helicopter rotor tip is in fact time 
dependent. It takes the form U = UR + V sin fit. This time dependence will be ignored in 
the derivation of Theodorsen’s lift and moment coefficients. 


Substitution of the Fourier expansions of 0, 37 , and Af into Equation (9) gives 


h r ( A1 '° + f ,/Sn*. = vo 0 

r + E ("». + ■»«*. (r x ‘) + d ~t) e '" a 


dh, 

dt 


(13) 


Equation (13) can be solved for to give 


**u(o»+h 4 k) 




Ar„ = 


27 TU (fl n +^»(j -»,.)+ ft^f ) 

J,°° \/lt^| c_,nnz " /, ' dx 0 


(14) 


Since AT n depends only up an n-th harmonic of 0 and -A, it is possible to write the 
the perturbation velocity potential as the linear combination of the velocity potentials for 
each harmonic. 



where 


^=u[o 0+ i- d -^)^A^ 


+ 2 irU 


/. + i*.\ ir 

\ Udl) r/Kii*. 


” — dz() 


, ,, f. tnfl (x \ „ 1 dh n r 

*• =u [®" + 7 T (5 ' «) + u rffj vTr * 

_ rr [. mfi / 1 \ ,, 1 dh„ 

+ I*V[>.+ — — jx 


fOO, _ 1 \/l-X 3 ^/*i*~* J 

Jl tan l^V- rf3r O 
f™ \f^* -' nUl, ' /U d * 0 


Since small perturbations were assumed, the lifting pressure distribution on the airfoil 
is linear in and in Therefore, the lifting pressure can be written as the linear 

combination of the lifting pressures calculated for each harmonic, 


*p=-*toZ(-ar + u -te ) 

n — fl x ' 


and the lift and moment coefficient for an arbitrary heaving and pitching motion can be 
written as the linear combination of Theodorsen’s lift and moment coefficients for each 
harmonic. 


c l ~ r -l, -* f/ - 

n. I 

V 

C‘rn = C (t j H ■+ ^ Cni„ 


The lift and moment coefficients used in this study were obtained from Reference 3 and 
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are as follows. 


and 


n/? 

V 


c l„ = c /„ ^0 + Jj + Co (1 - Cpa) 

«/. = C L - Co (3/4 - Cp a ) J O n + 


G d 2 h 


n + Fdh n 


nnu dt 2 U dt 


+ (c,(3/i- C; , a )- + — )(~4n/?„) 

c„ ( (do n i i d 2 h . i d 2 0„\l 

+ 4 \ \ dt + n/? ") U + V 2 dt 2 ColC, '“ 1/2 ^ (/ 2 dt 2 Jj 


'trit, 


-Cl., {c, a - 1 /4) (<?o + —£■ + c 0 (3/ 4 - Cp a ) ^ 


fid 


wi . . . n/? 0 l 

' f (C/.a “ 3/4) — j 


C m„ — c l„ 


(c,,„ - 1 /4) ( (f - c„ (3/4 - c,.„) ° jP- ) I. 4 ■ 


G d 2 h n 
nnu dt 2 


+ u t? * ( f ■ > 1 13/4 ‘ o " ;s ) ( ' T " ni> " ) ) 

+ J ((■>« - 3/4) -t- Wn) - + - 1/2) --f p 

F and G are. functions of the reduced frequency k, ■ 


k = 


Cj,nU 

2(7 


and are given by 


F + iG = 


/,° 


xu_ e -ifcx„ rfro 


r/x.i- 


f°° _£ii±i_, c -i*:xi, Jt 


( 19 ) 


( 20 ) 


( 21 ) 


( 22 ) 


d,(fc)-»r,(fc) 


J. (fc) + Vo(A-) - i (K, (k)-Jo(k)) 
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The J,'s and the K,*s are «-th order Bessel functions of the first and second kind. 

Plots of the lift and moment coefficients for sinusoidal pitching, sinusoidal heaving, 
and combined pitching and heaving motion are plotted in Figures 4 to 7. The experimental 
data was taken from Reference 4. These figures show that there is some discrepancy in 
the amplitude and phase of the lift and moment coefficients. This may be due to viscous 
effects which are ignored in Theodorscn’s theory. 


3.2. Equation of Motion for tho Frce-Tip 


The pitch angle measured with respect to the disk plane is given by the differential 
equation 


d 2 0 bddO kfi V* Moments 

j. _ ... 4 .S.Q ~ *=< 

dO Jdt J J 


(23) 


where 

bg ^mechanical damping coefficient 
k$ =mcchanical spring constant 


The moment forces in the pitch angle equation include 


9<i ^constant control moment 


A/, i ^aerodynamic moment at the quarter chord 

=c n „. ?c. 2 , 


M j =momcnt due to friction 




<i* 


A// e -feathering moment 

= - n 2 ("ijlj < sin 0 


M f =acrodynamic spring 
= - Ac ci g c 2 

M,ih “moment due to drag and dihedral 
=c,iqc 0 K ' n P o/.costf 
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The pitch angle measured with respect to the disk plane (Figure 3) is simply given in terms 
of the local angle of attack, or, and the inflow angle, <&, by 


6 = ft - <j> 


(24) 


The lift and moment coefficients given in Equations (19) and (20) indicate that they can 
each be broken into two terms, the first which is dependent only upon a, %f, and and 
the second which is dependent only upon <t > , , (i . and ■£. Letting 


c/ =cj(a) + c/(^, ~ , 0) 
c m =c m (o) 4 c m (<t>,~,0) 


(25) 


and substituting into Equation (23) gives 

d 2 a b» da k# c ( [a)qc 2 Ac c m (a)qc 2 

TP+lTt + 7 Q + — j T~ 

-f mj 2 ) q ( , d 2 d> be d$ k e J 

4- a - — t 4- — <p 

J J dl* J dl J 

ci(d>, gi»ff)< 7 C(.Ac c (<4, $ , /j)<?Cn Ac 

J + j 

n 2 (m/^ + m a f 2 ) A// A/,,/, 

f j ' s+ 'j ' T 

or defining new operators, L and At, 

£ (“) - — j — 1 ~r 


(26) 


(27) 

(28) 


where 

A At = M,ih + A 1j 

The function A At is not a linear function of a and its derivatives. It is, however, small 
with respect to the aerodynamic spring and can be ignored as a first approximation. 

By expanding a into a Fourier scries and noting that the operator £ is of the form 
4- B(t) -J* 4 C(l) a. Equation (27) can be written as 


Ou£(l) + ^ /l„£{cos»nt) 4 B n t(s\nnfU) 


(29) 


j 
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Both £(/(<)) a nd M(/( <)) can b® expanded into Fourier scries, 

30 3C 

£(1) =/!oo+ ^Om C03 + ]C ^omsinmflt 

me I m= I 

00 oo 

£(cosnnt) =/l c .o+ cos mill -t- £ B Chm sinmIlt 

m~l m^l 

00 oo 

£(sinnfR) =/l,„o + ^ A,„ m cosmnM ^ /i #i> „,sinmn/ 

m = J m ~ 1 

oo oo 

.M=.Mo + -Mem cos mill + M, m s\nrnHt 


m=l 


me 1 


( 30 ) 


Substitution of these Fourier expansions into Equation (29) and solving for the Fourier 
coefficients of a gives an infinite system of linear equations. 


(31) 


Two different computer programs were developed to solve for the Fourier coefficients 
of alpha. The first program, UNSTEADY was developed for a free-tip rotor in forward 
flight. The number of harmonic terms for o. A, d. and was limited to ten. This reduces 
the infinite system of linear equations for the Fourier coefficients of o to twenty-one. The 
coefficients were then solved by using Gauss elimination. Once a solution was determined, 
the right hand side of the equation of motion was then corrected by AM, and the harmonic 
coefficients were redetermined. This procedure was iterated until the solution converged. 

The second program was UNSTHOVR. This program determines the response of the 
free-tip to an air jet. The number of harmonics was limited to 30 which reduces the 
number of linear equations to sixty-one. The method of solution was identical to that of 
UNSTEADY. 


Mo 


ttoAijo + ^ An A Ch 0 + Xrf Bn A )u 0 = ~J- 


n = 1 

oo 


n= I 

DO 




TX~ I 

OO 


n I 

cc 


OoBum 4 ^ B, u „, 1 ^ - - j 


n = I 


»i - 1 
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4. Results 



4.1. Small Scale Model 

4.1.1. Rotor parameters and numerical model description 

A small scale model with a rectangular free-tip was studied in both hover and forward 
flight. The rotor parameters selected for the analysis are 

Radius = 8.40 ft 
Chord length =0.50077 ft 
Free-tip span = 0.84 ft 
Moment of inertia = 0.388 x ]0“ :< .slug ft 7 / ft 
Average tip speed = 700.0 ft /sec 
Forward flight speed =0.0 and 209.4 ft /fee 
9o = 5 lb ft/ ft 
H =0.0 and 0.3 
Airfoil : V23010-1.58 

The parameters which contribute to the feathering moment are 

mj =0.1133 x 10" 2 slug/ ft 
>»i« =0.1012 x 10" 3 slug/ ft 
//- 0.0075 ft 
l a =0.0840 ft 

Note that the parameters for the free-lip dimensions are two dimensional. The center 
of gravity of the free-tip is located at. its pitch axis. For a forward flight velocity of 
V = 209.4 ft /sec. the inflow and flapping angles were calculated by the Boeing Vertol 
B-G5 program. Two sets of data for rotor blade flapping and inflow angle distribution 
at the tip were generated bj B-G5 with CY/'rr ss ().()« and C\/a ss 0.007. Boeing Vertol 
B-65 is a three dimensional rotor performance prediction program for elastic blade rotors, 
employing the blade element theory (Reference 2). The first data set was computed with 
a conventional rectangular fixed tip rotor, and the second set with the free-tip, whose span 
length was 0.84 ft. When the inflow angle was computed with the B-65 for the frcc-tip 
rotor, it had an elastic inboard section, which is up to 90 % of the blade where the free-lip 
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it 
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starts. For fixed tip rotor, the entire blade was elastic. In both cases, the blade elasticity 
was prescribed in terms of the first elastic mode shape in flapwise bending and torsion. The 
free-tip and fixed tip have an identical planform. The inflow and flapping angles obtained 
from B-65 were then used in FTR.2 and UNSTEADY. These two sets of data are referred 
as fixed tip inflow set and free-tip inflow set respectively from now on. The resulting local 
angles of attack from the prescribed inflow angle and the rigid blade flapping, calculated 
by FTR2, B-G5, and UNSTEADY arc shown in Figure 8 for the fixed tip inflow set and in 
Figure 9 for the free-tip inflow set. 


It should be noted that the computer programs, FTI12 and UNSTEADY do not make 
iterations for the inflow angle. In other words, the resulting lift on the free-tip does not 
affect the inflow angle. This will limit the accuracy of the analysis. Even for the inflow 
angle, generated for the free-tip, the inflow angle contains inherent inaccuracy because the 
data was generated by B-65, which ignored the higher harmonic lift oscillations. 


For these two data sets, there is disagreement in the three solutions. This can be 
partially explained by the differences in the values used for the lift curve slope. This had 
its largest effect on the advancing side where the formula for the lift slope in FTR2 and 
UNSTEADY has a discontinuity in ^ y . However, this can not account for all of the 
discrepancies since both FT1I2 and UNSTEADY used the same lift curve slope. 

The additional differences may be explained by the differences in the lift and moment 
coefficients used in each program. In B-G5, Thcodorsen's unsteady lift and moment coef- 
ficients were used. The assumption was made that the constant and first harmonic terms 
dominated the Fourier expansions of the angle of attack, the heaving motion, and the 
flapping motion. No attempt was made to treat the harmonics separately, instead, they 
were all treated as if they were the first harmonic. 


In FTR2, the steady forms of the lift and moment coefficients were used. The damping 
coefficient was assumed to be the damping coefficient calculated for the first harmonic of 
UNSTEADY. Again, no attempt was made to treat the harmonic terms separately. 

In UNSTEADY. Theodorsen’s lift and moment are uniquely determined for each har- 
monic. 

To determine what effect these differences in angle of attack distributions might have 
upon the solution, consider the case of hover with an air jet. Since U is constant and 
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- 
dt 3 


the lift and moment coefficients take the form 
c/„ =C 0 ao 
c /« = CVi Of» + Z) n 


da n 


Ctriu ffo ®0 
Cm, = ffn a n "h ff n 


(ft 


tfftn 

(ft 


( 32 ) 


The operator Z acting upon the n-th harmonic gives 

£(«) = §+ ( J + (Ac On - K„) JC?,) ~ 


A'. , (AcC„ - lln)<ld 


(7 


)° 


or 


, d 2 o . b 0 da . A $ + A- u 

CW = 5Ti + — -rfi H 


(33) 


(34) 


where C„, £)„, W„. and /\„ are constants determined from Tlieodorsen’s lift and moment, 
and k a and f» t> are the aerodynamic spring constant and the aerodynamic damping, kg is 
a mechanical torsional spring and bp is a mechanical viscous damper. 

For HOVER, C„ and 11 „ in the above equation are replaced by 


For B-f>5, 


C n = Co 
II n=- 0 
D n = D\ 

Kn = I<1 

C n =C, 
D n =Dr 
n n = Hr 
K n = K x 


(35) 


for n - 1 , oc 


(36) 


for n = 2,oo 


If the mechanical damping and mechanical spring are ignored, the equation for the 
n-th harmonic of a is 

d 2 a„ 


dt : 


+ ~ ~~~ + -jo - |At„| sin (nflt + <f> n ) 


(37) 
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, 5- - ;rr, 



T 


Since U is constant for hover, the terms A, u m » A tu m, B,„ m , and B eK m in Equation (31) 
are equal to zero unless n " m. 

The amplitude of a n is given by 

| A n |= =JJ I = (38) 

1 1 v/(ifc-n 2 n 2 )2 + (6nn) 2 

If the forcing terms for HOVER, B-65, and IJNSTHOVR are equal, it is possible to deter- 
mine the ratio of the A n calculated by UNSTHOVR with respect to the A n calculated by 
HOVER or B-65. In UNSTHOVR the flapping motion is identically set equal to zero. The 
results are shown in Table 1. Compared to the coefficients calculated by UNSTHOVR, 
HOVER overestimated A„ for the first four harmonics; B-65 overestimated the second 
through fourth terms. However, for the fifth through tenth harmonics, both FTR2 and 
B-G5 badly underestimated A n . 

The forcing function M is also affected by the treatment of the higher harmonic 
terms. The function M includes lift and moment terms which are due to the inflow angle, 
the heaving motion, and the flapping angle. The lift and moment due the steady portion 
of the inflow angle cancel out the lift and moment due the steady portion of the heaving. 
However, this is not true of the higher harmonic terms. 

In HOVER no correction was made to >1 to allow for these higher harmonic lift and 
moment terms. B-65 treated these corrections as. if they were all first harmonic. For no 
mechanical damping and no mechanical spring. .M becomes 

d 7 o h d<? k n 2 (»» jlj - #h„/, 2 ) g„ AM 

d7 2 " Jdt* 1 ' ‘ .7 J + ~J~ 

where 6 and k are the damping and spring coefficients derived from the c/(d>,/3, -j|) and 
f; ). Table 2 shows the ratios of the forcing function amplitudes and the harmonic 
coefficients for the angle of attack, a. Note that there is a sizable correction to the first 
harmonic term of M in UNSTHOVR. 

4.1.2. Free-tip response to an air jet during hover 

The free-tip response to a vertical air jet during hover is now investigated. The air 
jet is modeled as a change in inflow angle. For simplicity, inflow angle change across the 
air jet was assumed to be a half sine wave. Note that the frequency of the excitation due 
to the air jet can be changed by changing cither the rotor angular velocity or the width 
of the air jet. In this study, the width of the air jet was selected to be a variable so that 
all rotor parameters can remain constant. Figure 10 shows the azimuthal lift distributions 
computed by HOVER, and Figure 11 by UNSTHOVR. The results from UNSTHOVR also 
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include the lift distributions when two different levels of damping were added. For these 
two cases, HOVER and UNSTHOVR, the inflow angle was prescribed as follows. 


^ =0° for 0° < rjj < 90° 

4> =3 0 sinCnt for 90° < tl> < 120° 
<£ = 0° for 120° < 0 < 360° 


There are significant differences in the solutions. The peak lift value for the HOVER cal- 
culation is much higher than UNSTHOVR. This occurs because the aerodynamic damping 
coefficient calculated for the first harmonic is less than one half of the aerodynamic damp- 
ing coefficient calculated for the higher harmonics as shown in Table 1. Since the program 
HOVER ignores the higher harmonic terms, the lift distribution behaves as a simple un- 
derdamped oscillator. This is not the case with UNSTHOVR where each harmonic term 
has its own aerodynamic damping and spring coefficients. Note especially how the fifth or 
sixth harmonics are not damped out in Figure 11. 

For Figures 12 and 13 the inflow angle was changed as follows. 

<f> - 0° forO" < 0 < 90" 

4> — 3° sin 12Qt for 90 n < t!* < 1 05 
4> - 0° for 105” < 0 < 3G0" 

Again, the result of UNSTHOVR includes the additional lift distributions with two levels 
of damping. In spite of the high excitation frequency, 12ft. Figure 13 shows that the 
relatively low harmonic oscillation is not damped out. namely, the fifth harmonic. The 
lift distributions in Figures 11 and 13 suggest that the natural frequency of the system 
is between 5ft and 6fi. From Tabic 1, it is also possible to determine that the natural 
frequency for this system is somewhere between 5ft and 6ft. Figure 1-4 shows a difference 
between the square of the natural frequency. of the system, -f , and harmonics of the 
rotor angular velocity to the square, (nft) 2 . Interpolation of the data suggests that the 
natural frequency of the system, is approximately f>.2fi - 433.33 rad/sec. For both cases 
increasing the mechanical damping also decreases the amplitude of lift oscillations. 
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4.1.3. Free-tip response in forward flight 

For forward flight, the effects of friction, mechanical damping and mechanical spring 
on the azimuthal lift distributions are investigated. The computations were carried out 
with two different inflow sets described previously. The fixed tip inflow set was used for 
Figures 15 to 17, and the free-tip inflow set for Figures 18 to 20. In all cases, there is a 
noticeable lift peak on the retreating side around t/> = 300° to 330°. This peak is caused by 
the dynamic pressure build-up. Since the dynamic pressure on the retreating side is small, 
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the aerodynamic spring on the free-tip, which determines the capability of the free-tip 
response, is also low. Therefore, the free-tip fails to respond to the fast dynamic pressure 
build-up beyond tp = 270 3 , which results in the lift peak. This is the g-effcct discussed 
in Reference 1 and 2. For both sets of inflow angle and flapping angle data there is a 
large amount of variation in the lift force on the advancing side. This variation can be 
smoothed out by decreasing the coefficient of friction (Figures 15 and 18) or by increasing 
the mechanical damping (Figures 16 and 19). Increasing the mechanical damping has the 
undesirable efTect of making the free-tip slow to react and causes the peak to peak lift to 
increase. 

Figures 17 and 20 indicate that increasing the mechanical spring constant helps to 
decrease the peak to peak lift with the fixed tip inflow set, since the ( 7 -effect is the dominant 
in the lift distribution and adding a mechanical spring increases the natural frequency and 
enables the tip to respond faster. With the inflow angle distributions generated by B- 
65 with the free-tip, the (/-effect does not seem to be significant because of its inflow 
pattern, and adding mechanical spring causes larger oscillations on the advancing side due 
to underdamping, which was predicted by the unsteady lifting surface theory (Reference 
1 and 2). There is no appreciable smoothing of the higher harmonic terms of the lift 
distribution by the inclusion of mechanical damping. 

Although the rotor parameters were identical for the two small scale cases, the az- 
imuthal lift distributions calculated for the two sets of inflow and flapping angles are 
dissimilar due to the difference in inflow pattern. The mean lift level calculated with the 
free-tip inflow set is higher than that calculated with the fixed tip data. The lift distribu- 
tion with free-tip inflow set has more higher harmonic components than the lift distribution 
with the fixed tip inflow set, although its peak to peak lift is smaller than the peak to peak 
lift calculated for the fixed tip inflow set. 

In Figures 21 and 22, the pitch angle of the free-tip computed with the fixed tip inflow 
set and free-tip inflow set arc presented. The fixed tip inflow set produced a much larger 
peak value for the pitch angle. For the fixed lip inflow set., vortices shed from the tip of the 
inboard section do not exist. In the free-tip data, the effect of these vortices is included. 
The inclusion of these shed vortices causes an upwash component to be added into the 
inflow angle at the free-tip. For this model the upwash component was not large enough 
to make the inflow angle positive. It did. however, reduce the magnitude of the inflow 
angle on the retreating side. This caused the pilch angle to remain positive and to have a 
smaller peak value. 
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4.2. Full Scale Model 
4.2.1. Rotor parameters 

A parametric study was performed for a full scale model with a swept free-tip. The 
parameters of interest were the sweep, the constant control moment, a constant moment 
coefficient at 25% chord, the mechanical spring constant, mechanical damping, and friction. 
The criteria for success was a constant lift distribution. 

The rotor parameters used in this study are 

Radius = 28.0 ft 
Chord length = 1.2657 ft 
Free-tip span = 2.8 ft /sec 
Moment of inertia = 0.02 slug ft 2 /ft 
Average tip speed = 650 ft /sec 
Forward flight speed -- 325 and 95 ft/sec 
9u =40 lb ft/ ft 
H — 0.5 and 0.15 
Airfoil : V23010-1.58 


The parameters for the feathering moment are 


m f =0.1649 x 10 
m n =0.4415 x 10 
/ n = 0.28 ft 
t f =0.075 ft 


l 

-2 


slug' ft 
slug! ft 


Note that the parameters for the free-tip dimensions are two dimensional. The center of 
gravity of the free-tip was again located at its own pitch axis. Two sets of data were used 
for the flapping and inflow angles. For the first, set of data, the forward velocity was 325 
ft /sec with the advance ratio of (i = 0.5. The inflow angle was calculated by B-65 for a 
fixed tip with C't!o 0.08 and Cx/o ss 0.007. No flapping angles were available. 

For the second set of data, the forward flight velocity was 95 ft /sec. with the advance 
ratio beingp = 0.15. The Happing and pitching motions were obtained from flight test data 
for a S icons ky S-58 helicopter with a fixed tip rotor (Uefcrence 5). The thrust coefficient 
based on the weight of the helicopter, 11,500 lb was Ct/° 0.08. The inflow angle was 
calculated from the experimental lift and the pilch angle. 
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4 . 2 . 2 . Free-tip response in high speed forward flight 

The lift distributions calculated with the first set of data with n = 0.5 are presented 
in Figures 23 to 29. The lift distribution shows a large amount of oscillation on the 
advancing side. There is also a large peak in the lift distribution on the retreating side. 
This is the 9-effect discussed previously. Increasing the friction coefficient to 0.2 x 10~ 2 
reduces the amplitude of the lift distribution (Figure 23). This behavior is similar to what 
was observed in small scale model. Since the friction is dependent upon the sign of 
it acts as a damping term. When the motion is underdamped, additional friction helps 
to smooth out the lift distribution. When the motion is overdamped, additional friction 
contributes to the oscillatory lift distribution by further suppressing the capability of the 
free-tip to respond to the perturbations. 

The inclusion of mechanical damping smooths out some of the higher frequency oscil- 
lations (Figure 24). As was the case with the small scale model, large mechanical damping 
prevents the free-tip from reacting quickly and causes the peak to peak lift to increase. 

Inclusion of a mechanical spring has little or no effect on smoothing the higher fre- 
quency oscillations on the advancing side (Figure 25). However, it did have a major effect 
on flattening the lift distribution on the retreating side, since the mechanical spring helps 
to minimize the 9-effect where the aerodynamic spring term is low. 

Decreasing the effective pilch axis offset from the aerodynamic center of the free-tip 
had a drastic effect upon the lift distribution. Decreasing the effective pilch axis decreases 
the aerodynamic spring and damping. This allows the free-tip to oscillate with a larger 
amplitude as indicated in Figure 26. 

The major effect of lowering the constant control moment was to lower the average 
lift ( Figure 27). There are some small changes in the shape of the lift distribution, but 
very little smoothing occurs. 

The inclusion of a constant aerodynamic moment coefficient at the 25% chord pre- 
dominantly changes the lift on the advancing side (Figure 28). The positive c m raises the 
lift and the negative c,„ lowers the lift. There is no smoothing effect. 

It appears that to obtain a constant lift distribution, both mechanical damping and 
a mechanical spring must be added. The mechanical damping is used to smooth the 
lift distribution on the advancing side; the mechanical spring is used to smooth the lift 
distribution on the retreating side. Both the mechanical spring and damping must be 
added since the mechanical damping causes a larger 9-effect on the retreating side and 
the mechanical spring causes more lift oscillation :>n the advancing side. The combination 
of mechanical damping and mechanical spring are shown in Figure 29. This combination 
does indeed smooth the entire lift distribution and decrease the peak to peak lift, since it 
increases the natural frequency of the tip, without changing the damping characteristics. 
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4.2.3. Free-tip response in low speed forward flight 


The lift distributions calculated with the second set of data with fi = 0.15 are plotted 
in Figure 30, with three different levels of mechanical damping. The experimental set of 
data produced a much flatter distribution than the previous set of data for the high speed 
flight for all three levels of damping. This occurred because the higher order harmonic 
terms in the inflow angle, the downwash velocity, and the forward velocity are smaller than 
those in high speed data set with /j. = 0.5. 
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5. Conclusions 




Two computer programs, UNSTEADY and UNSTHOVR were developed for the inves- 
tigation of the free-tip rotor response duri.ig forward flight and hover. They both incorpo- 
rate Theodorsen’s lilt deficiency functions, uniquely defined for each harmonic component 
of the free-tip pitching and heaving motion. 

A large pitch axis offset of the free-tip from its aerodynamic center increases the 
aerodynamic spring, which improves the free-tip response. The addition of the mechanical 
spring reduces the (/-effect on the retreating side. Mechanical damping must be added 
to minimize the oscillation on the advancing side. Neither the mechanical spring nor 
the damping should be large since the mechanical spring increases the amplitude of the 
oscillation on the advancing side and mechanical damping increases the ^-effect. This result 
agrees with those cited in Reference 1 and 2. The calculations indicate that to achieve 
a constant lift distribution, both mechanical damping and a mechanical spring must be 
added. 

Even with the higher harmonic terms on the advancing side, the lift distribution for 
the free-tip is much more constant than the lift distribution for the conventional fixed tip. 
It does not have a large negative lift on the advancing side. This alone should reduce some 
of the oscillation associated with the fixed tip rotor. 

The higher harmonic terms of the angle of attack, the heaving motion and the forcing 
function in Theodorsen’s theory were ignored in FTR2 and B-65. These progarms treated 
the lift and moment as if there was only a first harmonic motion. 

It should also be noted that the inflow and flapping angles played a major role in the 
azimuthal variation of the lift coefficient as indicated by the differences in the azimuthal 
lift distributions of the two small scale models. Both models shared the same physical 
dimensions and the same forward velocity, however, the inflow and flapping angles were 
different. The lift distribution calculated with the free-tip inflow set had a higher average 
lift, a smaller peak to peak lift., and larger higher harmonic terms than the lift distribution 
calculated with the fixed lip inflow set. 

The free-tip calculation of the inflow and flapping angles by B-G5 treated all harmonic 
terms as if they were a first harmonic. In UNSTEADY and UNSTHOVR, the computer 
program was not iterated for the resulting lift and inflow angle, therefore, it is impossible 
to predict how the correct treatment of each harmonic term would affect the inflow and 
flapping angles and, hence, the azimuthal lift distribution. This should be seriously con- 
sidered when examining the full scale results. The full scale inflow and flapping angles are 
either taken from a program which had a fixed tip or they are taken from experimental 
data fo. a fixed tip rotor. The wake relaxation of the total rotor system with the correct 
treatment of higher harmonics should be performed to solve this problem. 
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Table 1. Ratios of the harmonic coefficients 
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COEFFICIENT RATIOS 

HARMONIC 

*r 

-V - (" n ) J 

M'*'l 

ld!!U 

M‘ a M 

1 

54.07 

226,947 

0.940 

1.000 

2 

135.58 

192,852 

0.880 

0.941 

3 

180.87 

146,599 

0.823 

0.892 

4 

210.18 

87,830 

0.812 

0.905 

5 

230.58 

16,128 

1.248 

1.514 

6 

245.58 

-68,856 

5.195 

4.471 

‘ 7 

256.59 

-167,395 

2.310 

2.024 ~] 

j 8 

265.20 

-279,695 

K659 

T549 1 

! 9 

r 272.97 

-405.915 

1.434 

1.372 j 

10 

277.39 

-546,176 

1.323 

1.282 1 


A is the value calculated by HOVER 
A is the value calculated by 13-65 
A^ is the value calculated by UNSTHOVR 
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Figure 1. Free-tip rotor schematics 
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Figure 16 . Effect of mechanical damping, fixed tip inflow set 
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Figure 21. Pitch angle of the free-tip with fixed tip inflow set 
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Figure 22. Pitch angle of the free-tip with free-tip inflow set 
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Figure 28. Effect of the aerodynamic moment coefficient 
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Appendix A 


. Input Data for UNSTEADY 


Sample Data Set 


$SUSIE NMAX-20. EPS-3. 0000 1 

N-144.J-. 02. DELC-0. 232. CPA-3. 125, C-1.2657, 

QUSO-1 436.36, VSOUND— 1 143.75. V-192.3. OR-630.0. RADIUS-28. 0. 

UASS-0. 18634. QZERO-4O.0. A1-0.O, B1-0.O. CM0-O.O, PHIZ-7. 43 125. 
BETA0-O.0.D1— 5. 63269, 02-3. 75632. D3-0. 94621. 04— 0.52603.05— 0.46910, 
D6-0. 09250, 07-0.16116, DO- 0.14434, 09—0.05079. 010—0.16852, 
FI-0.62908, F2— 0.53383, F3- 0.32642. F4- 0.36417, F5-0. 14365. 

F6— 0.11166, F7- 0.09747. F8-0. 20000. F9— 0.00473. F10— 0.CO7CO. 
BS1-.0. BS2-0.0, BS3-0.0, BS4-0.O, BS5-0.0, 

BS6-0.0. BS7-0.0. BS8-0.0, BS9-0.0. BS1O-O.0. 

BC1-0.0. BC2-0.O. BC3-0.0, BC4— 0.0. BC5- 0.0, 

BC 6-0.0. BC7-0.0, BC8-0.0. BC9-0.0, BC10-0.0, 

CLAO-3.552, DIVMAC-0.75, 

CO0-0. 0034. C01— 0.0102. C02-0.384. VWCRIT-0.83. ACRIT-33.. 
DCDDW-.05, MUL-0.001. BOOTH-O.O, K THETA-0.0, 

MA-0. 004415, UF-0. 016490. LA-0.28. LF-0.075. 

AP0- 22.71, API— 2.23, BP1-13.28. TIT— 8.0. TL0C-25.20. 

TIPL-2.8, CAWMA-0.0. THPRE-0.0 $EH0 


NMAX 

EPS 

N 

J 

DELC 

CPA 

C 

QMSQ 

VS0U1ID 

V 

OR 


Input Data Description 

maximum number of iterations 

limit on a n+ i - a n (rad) 

number of computation points 

moment of inertia of Free-Tip ( slug ft/ ft) 

effective pitch axis offset, Ac 

location of the. pitch axis measured from the leading edge 

Free-Tip average chord length (ft) 

dynamic pressure/(Mach number) 2 (lb/ ft 2 ) 

sonic velocity (ft/sec) 

helicopter forward flight velocity (ktu) 

average tip speed (ft /sec) 



RADIUS 
MASS 
qZERO 
Al. B1 
CMO 

PHIZ. Dl-DIO 
F1-F10 

BSH1-BSN10. 

BCS1-BCS10 

CLAO 

DIVMAC 

CDO.CD1.CD2 

VMCRIT.ACRIT 

DCDDM 
MUL 
BDOTH 
KTHETA 
MA. MF 

LA. LF 

APO. API. AP2 
TW 


rotor radius (ft) 
mass of the Free-Tip (slug/ ft) 

constant control moment (ft lb/ ft) 

cyclic control moment (ft lb/ ft) 
constant moment coefficient 

harmonic coefficients which define the inflow angle distribution 
(deg) 

harmonic coefficients which define the flapping motion (deg) 


lift curve slope of the tip section at = 0 

lift divergence Mach number (See listing.) 

coefficients which define the drag coefficient as a function of a 

critical a and Mach number used to define the compressibility 
drag increment 

at critical Mach number, VMCRIT 
friction coefficient 

mechanical damping factor (/6 sec 2 / ft) 
mechanical spring rate (lb ft/ ft) 

forward and aft masses which contribute to the feathering mo- 
ment (slug/ ft) 

forward and aft moment arms which contribute to the feather- 
ing moment (ft) 

collective and cyclic pitch control of the inboard section, meas- 
ured at the blade root in the shaft axis system (deg) 

blade twist, negative down (deg) 
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JtaJWT Mli 



/ 

/ . 



TLOC 

TIPL 

GAMMA 

THPRE 


/ 


TT . 
/ i 






si' ■ \/ 


radial location of the Free-Tip c. g. (ft) 
Free-Tip span length (ft) 

Free-Tip dihedral angle, positive down (deg) 
mechanical spring pretwist angle (deg) 
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Appendix B. Input Data for UNSTHOVR 

Sample Data Set 


<- 



tSUSIE NMAX-40. EPS-0. 00001 
N- 144. J-. 000388. 0ELC-.1 3. CPA-. 125, 0.56077 
QMSQ-1439.37, VSOUHD— 1 137 . 73, V-0.0. OR-700.0, RADIUS-8.4. 
UASS-0. 01157, OZERO-0 .0. A1-O.0, Bl-0.0, CW0-O. , 

CLA0-3 . 552 . DIVMAC-0. 75, 

CDO— 0.0101 . CD1— 0.008, CD2-0.075, WCR IT-0. 706. ACR1T-32.2. 
DCDOM-0. 105. WUL-0 . 001 , BDOTM-0. 0, KTHETA-0.0. 

MA-0. 0001012, MF-0.001133, LA-0.084, LF-0.0075. 

AP0- -6.48. API-0.0. BP1-O.0, TW—9.450, TL0C- 7.56, 

TIPL-0.84, CALO4A-0.0. THPRE-0.0, PHI A-3 .0. PH IB-0. 0 , PST ART-90 . . 

P END-1 20.0 SEND 


l 




UMAX 

EPS 

M 

J 

DELC 

CPA 

C 

QMSQ 

VS0U1JD 

V 

OR 


Input Data Description 

maximum number of iterations 

limit on o n .| a„ (rad) 

number of computation points 

moment of inertia of Free-Tip (.si tig ft, ft) 

effective pitch axis offset, Ac 

location of the pitch axis measured from the leading edge 
Free-Tip average chord length (ft) 
dynamic pressure/(Mach number) 2 (lb/ ft 2 ) 

sonic velocity (ft /nee) 

helicopter forward flight velocity ( kts ), must be 0 
tip speed (ft /see) ' 

01 
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, • >-* / ‘ 

/ . //} ■; ■ / ■' . / 

\ 

“ ■' ° . 0 V ; 


/ : 

/ . 

RADIUS 

rotor radius (ft) 

. , 

' rp • 

MASS 

mass of the Free-Tip (slug/ ft) 

■'-1 

* ••/ 

qZERO 

constant control moment (ft lb/ ft) 

\:hL.. 

A1.B1 

cyclic control moment (ft lb/ ft) 


CMO 

constant moment coefficient 


CLAO 

lift curve slope of the tip section at M = 0 

■ / 

DIVMAC 

lift divergence Math number 


CDO.CDl ,CD2 

coefficients which define the drag coefficient as a function of a 


VMCRIT.ACRIT 

critical a and Mach number used to define the compressibility 
drag increment 


DCDDM 

at critical Mach number, VMCRIT 

; 

MUL 

friction coefficient 


BDOTH 

mechanical damping factor (lb sec 2 ’ft) 

/' 

KTHETA 

mechanical spring rate (lb ft/ ft) 


MA.MF 

forward and aft masses which contribute to the feathering mo- 
ment (slug/ ft) 


LA. LF 

forward and aft moment arms which contribute to the fcather- 
- .. ing moment, .(ft) . . ■ 

■ • f 

APO.AP1.AP2 

collective and cyclic pitch control of the inboard section, meas- 
ured at the blade root in the shaft axis system (deg) 

/ 

TW 

blade twist, negative down (deg) 

' / 

TLOC 

radial location of the Free-Tip c. g. (ft) 

f 

i : 

TIPL 

Free-Tip span length (ft) 

-'A 
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GAMMA 

Frcc-Tip dihedral angle, positive down (deg) 


THPRE 

mechanical spring pretwist angle (deg) 


PHIA 

amplitude of half sine inflow angle (deg) 


PHIB 

constant inflow angle (deg) 


PSTART 

azimuthal angle at which the half sine inflow angle begins 

(deg) 

PE11D 

azimuthal angle at which the half sine inflow angle ends 

(deg) 
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